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In this paper we consider the effect of surface heterogeneity on the slippage of fluid, using two 
complementary approaches. First, MD simulations of a corrugated hydrophobic surface have been 
performed. A dewetting transition, leading to a super-hydrophobic state, is observed for pressure 
below a "capillary" pressure. Conversely a very large slippage of the fluid on this composite interface 
is found in this superhydrophobic state. Second, we propose a macroscopic estimate of the effective 
slip length on the basis of continuum hydrodynamics, in order to rationalize the previous MD 
results. This calculation allows to estimate the effect of a heterogeneous slip length pattern on 
the composite interface. Comparison between the two approaches are in good agreement at low 
pressure, but highlights the role of the exact shape of the liquid-vapor interface at higher pressure. 
These results confirm that small variations in the roughness of a surface can lead to huge differences 
in the slip effect. On the basis of these results, we propose some guidelines to design highly slippery 
surfaces, motivated by potential applications in microfluidics. 

PACS numbers: 83.80.Fg, 68.35. Np, 68.10c 



I. INTRODUCTION 

The question of the hydrodynamic boundary condi- 
tion (h.b.c.) for simple fluids has received much atten- 
tion recently, and a number of experimental studies have 
addressed this problem QSHSHSHSSlIi 
IH 13 El 111 EH Usually, when a simple liquid flows 
over a motionless solid wall, its velocity is assumed to 
be equal to zero near the wall. This assumption is a 
postulate of hydrodynamics. Although it is widely ver- 
ified for macroscopic flows of simple fluids, its validity 
at small scales remains an open question. Wall slippage 
is usually described in terms of an extrapolation length, 
usually denoted as the slip length, b : this corresponds to 
the distance between the wall and the position at which 
the linear extrapolation of the velocity profile vanishes. 
Analytical calculations and numerical simulations have 
shown that two parameters mainly affect the h.b.c: the 
roughness of the surface ^3 an( A the solid-liquid inter- 
action 0, 03> 0| . While roughness is usually expected 
to decrease slippage, hydrophobocity of the solid surface 
has been found to enhance the slip length. ^From an 
experimental point of view, many studies have been car- 
ried out recently but they do not allow to draw a simple 
and unique conclusion concerning the effect of roughness 
and solid-liquid interaction on slippage, as a huge dis- 
persion in the results is observed for seemingly similar 
systems. For instance, on "smooth" hydrophobic sur- 
faces, slip lengths b were obtained that range from a few 
nanometers PJ up to a few micrometers Even at a 
qualitative level, some results appear to contradict each 
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other. To cite an intriguing point as an exemple, rough- 
ness has be found either to increase Q or to decrease 
the slip length ^2- ^ is therefore important to clarify 
this question from a fundamental point of view, in order 
to understand the origin of this important variability of 
results. ^From a more "practical" perspective, the pos- 
sibility for a liquid to slip over a solid surface can also 
have important consequences for flows in microchannels. 
A large slip length reduces dissipation at the surfaces and 
therefore affects the pressure drop in a given channel geo- 
metry. It also may help in reducing the hydrodynamic 
dispersion due to the velocity gradient between the center 
and the channel boundary, with important applications 
for exemple in chemical or biological techniques of analy- 
sis (micro-chromatography). Finally, slippage is also an 
important phenomenon for flows in mesoporous media, 
where the understanding of the effect of mechanical dis- 
sipation is essential for applications such as dampers [2(| . 

In previous work |2l| , we have shown by Molecular Dy- 
namics (MD) simulations that the slip length is strongly 
affected by the roughness of the substrate for non wet- 
ting solid-fluid interactions, a point which could explain 
some of the experimental discrepancies. In particular, 
we have found, for a non wetting patterned surface, a 
thermodynamic transition between two different surface 
states as a function of the pressure : a "normal" state 
where the liquid occupies all the available volume, and 
a "super-hydrophobic" state, where the liquid flows over 
the roughness. In the "normal" state, the h.b.c. is well 
described by a no slip condition (the roughness decreases 
the slip length), whereas the "super-hydrophobic" state 
enhances a slip h.b.c, with slip lengths greater than in 
the smooth case (the roughness increases the slip length) . 

Our aim in this paper is twofold. First, using MD si- 
mulations, we present a more detailed study of the transi- 
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tion between the "normal" and the "super-hydrophobic" 
state, as a function of the roughness parameters, geo- 
metry and pressure. Note that the geometry studied 
here consists of grooves with different widths and heights. 
This provides an ideal geometry, from which one could 
expect to extract simple rules. Slippage over the sur- 
face is conversely studied in the different thermodynamic 
states. Second, we try to rationalize the effect of the 
surface heterogeneity on the measured slip length, as ob- 
served in the previous MD simulations. To this end, we 
develop a macroscopic description of the hydrodynamic 
flow on a microscopically heterogeneous surface, charac- 
terized by a spatially varying h.b.c. Beyond the compari- 
son to the MD results, this analysis allows us to propose 
an optimized surface heterogeneity, designed to obtain 
large slippage at the interface. 



II. MOLECULAR DYNAMICS APPROACH 

A. Phase properties of heterogeneous surfaces 

We consider in this study a liquid confined between 
two parallel solid walls. The bottom wall is decorated 
with a periodic array of grooves of height h and width 
a. Note that this differs from the system considered in 
reference plj . where a periodic array of square dots had 
been considered. Periodic boundary conditions are used 
in the directions parallel to the walls, the length of a cell 
of simulation is L. In our simulations, all the interactions 
are of the Lennard-Jones type : 
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The fluid and the solid atoms have the same molecu- 
lar diameter a and interaction energies e. The variation 



the 



(the index i,j — F,S refers to the fluid or 



solid phase) is a convenient control parameter that can 
be varied to adjust the surface tensions. The solid sub- 
strate is described by atoms fixed on a the (100) plane 
of an FCC lattice. We have worked with two values of 
cfs, cfs = 0.5 or cfs = 0.8, which correspond to con- 
tact angles (deduced from Young's law) of 9 — 137° or 
9 = 110°, respectively 01 • The simulations are carried 
out at constant temperature ksT / 'e=l. In flow experi- 
ments, the velocity component in the direction orthogo- 
nal to the flow was thermostatted ^tJ; m order to avoid 
viscous heating within the fluid slab. 

Figure ^ displays the evolution of the pressure with 
the distance d between the walls of the cell. The normal 
pressure is obtained from the average force on the sub- 
strates along the z direction. Since the number of atoms 
is kept constant, this plot is similar to a P — V isotherm 
(pressure versus volume) in the bulk. We have verified 
that changing the pressure by modifying the number of 
particules in the cell at fixed volume V leads to equiva- 
lent results. As for the case of dots |2J, this isotherm 



exhibits two branches, separated by a "Van der Waals 
loop" , indicating a phase transition between two possi- 
ble situations. As shown in figure where transverse 
views of the atomic configurations for different pressures 
are represented, this transition corresponds to a partial 
dewetting of the space between the grooves. At high 
pressures, the liquid occupies all the cell (see figure QJ, 
including the grooves. This state will be described in 
the following as the "normal" state. At lower pressures, 
partial dewetting is observed and a composite interface 
is formed. In the latter case, the space between the 
grooves is essentially free from liquid atoms (see figure 
0. This state will be therefore denoted as the "super- 
hydrophobic" state. Note that the range of distances for 
which the pressure decreases with d is thermodynami- 
cally unstable, its observation being merely an artefact 
of small size simulations. 

Qualitatively, the maximum of the pressure on the pre- 
vious Van der Waals loop corresponds to the capillary 
pressure P cap , for which the liquid starts invading the 
grooves. On the other hand, the minimum corresponds 
to the pressure at which the liquid reaches the bottom of 
the groove and adapts its contact angle in order to satisfy 
Laplace's law. 

A qualitative interpretation of this dewetting transi- 
tion can be given in terms of macroscopic capillarity. If 
we consider a system at fixed normal pressure Pn, the 
difference in Gibbs free energies between the normal state 
(case ^4) and the formation of a superhydrophobic state 
(case B) can be written as: 

Gb-Ga = L{(L-a+2h) (j LV cos9) + (L-a) -f LV +P N (L-a)h} 
The composite interface is therefore favored when 
P,r < P - P + -1lv-{cos9 + 1) 

1 N 1 composite — 1 cap ~r . \o j 

where we have introduced the capillary pressure P cap , 
defined as 



Pn. 



2jlvcos9 
L — a 



(4) 



Although a quantitative agreement can hardly be ex- 
pected in view of the small sizes in our system, we have 
checked that equation correctly predicts the general 
trends observed in our simulation, when the height or 
width or the grooves are varied. As an example, we have 
plotted the measured capillary pressure on Fig. |2J here 
defined in the simulations as the maximum pressure for 
which the dewetted state is stable (corresponding to the 
local maximum in figure . This plot corresponds to 
the limit in which the grooves have a large depth h, so 
that the "composite" pressure, defined in equation © 
should reduce to the capillary pressure, equation (J2J. As 
shown in this figure, the linear dependence of the pres- 
sure as a function of 1/(L — a) does confirm the validity 
of the macroscopic argument above. Moreover, the slope 
of the linear fit agrees within 15% with the predicted one, 
2j LV cos9 



3 



B. Slippage over heterogeneous surfaces 

We now turn to the main point of the present paper, 
namely the slippage of the fluid on the previous hetero- 
geneous surface. 

We have conducted parallel Couette flow simulations. 
A shear flow is induced by moving the upper wall with ve- 
locity U and the lower wall with a velocity —U. Typically 
we chose U = 0.3 in reduced Lennard- Jones units. We 
have verified that, in spite of the large value of the velo- 
city, the response is linear and the measured slip lengths 
are independent of the velocity U. Temperature is kept 
constant by thermostatting through velocity rescaling in 
the direction perpendicular to the flow. For flat walls 
|l8|. the slip length b is defined as the distance between 
the wall position and the depth at which the extrapo- 
lated velocity profile reaches the nominal wall velocity, 
v = U. In the presence of grooves, the same definition is 
used. Obviously the presence of grooves makes the choice 
of the wall position somewhat arbitrary. We choose to 
define the wall as the position of the bottom layer of the 
substrate. This choice will be seen to have only a negligi- 
ble influence on the results (since, as discussed below, the 
measured slip lengths are found to be much larger than 
the groove depth in the interesting "superhydrophobic" 
state). 

Typical measured velocity profiles are plotted in figure 
|3J As shown on this figure, the slippage strongly depends 
on the surface state : globally, while in the "normal" 
state the slip length is found to be very small (roughness 
reduces slip), in the "superhydrophobic" state the slip 
length is strongly enhanced compared to the smooth sur- 
faces case (with the same solid- fluid interactions). One 
may note that for the smooth top surface the slip length 
remains almost constant with the pressure, and its evo- 
lution is not affected by the transition between the "nor- 
mal" and the "super-hydrophobic" state. This observa- 
tion shows that for the confinement we consider here, the 
dynamics in the vicinity of the top wall is not affected by 
the bottom wall. 

We present in figure 0] the evolution of the slip lengths 
on the smooth wall and on the rough wall, with the 
distance between the walls, together with the measured 
pressure (here the number of molecules per cell is fixed) . 

The results in figure 0] were obtained for a velocity pa- 
rallel to the grooves. We have performed the same anal- 
ysis for a velocity perpendicular to the grooves and ob- 
tained very similar results, which typically differ by a 
numerical constant (the slip in the direction perpendic- 
ular to the grooves being less pronounced than in the 
direction parallel to the grooves, as expected). In order 
to simplify the discussion, we will only present in the 
following the results for a liquid flowing parallel to the 
grooves. 

As mentioned above, in the normal state, the rough- 
ness is found to reduce the slip. The slip length in 
the "normal" state is lower on the rough wall than on 
the smooth wall. On the other hand, in the "super- 



hydrophobic" state, roughness considerably enhances 
the slip length. In the "super-hydrophobic" state, one 
observes a strong dependence of the slip length with the 
pressure, an effect which can be correlated to the strong 
dependence of the shape of the liquid/ vapor interface 
with the pressure. Indeed, as the pressure decreases, 
the interface changes from an incurvated shape in the 
grooves to a flat shape laying on the top of the roughness. 

In order to compare with experiments, it is more con- 
venient to plot the slip length as a function of pressure 
(figure EJ. In this figure the pressure is made dimension- 
less using the capillary pressure defined in equation J3J. 
^From this kind of representation one may expect to be 
able to estimate the slip length for a specific roughness 
geometry, knowing the pressure of the system. From this 
figure, it is already obvious that small changes in the 
pressure induce very large difference in the measured slip 
length. 

We now turn to the influence of the geometry of the 
roughness on the slip length. We mainly focus on the role 
of the groove's width. Figure represents the evolution 
of the slip length as a function of the pressure for various 
widths of the grooves (keeping their height constant in 
this case). For a given pressure the slip length increases 
with the width of the grooves. This point can be ea- 
sily understood. A wider groove implies a larger area 
of liquid/vapor interface, leading to a lower friction at 
the interface. From this figure we also note that small 
variations in the geometry of the grooves can lead to 
very large variations in the slip lengths. This fact could 
explain the great variability observed in the experimental 
results. Small differences in the state of the surfaces used 
in experiments could lead to huge differences on the slip 
lengths. 



C. Conclusions and hydrodynamic description 

At this stage we can draw a few general conclusions. 
On the basis of MD simulations we have shown that the 
surface roughness is a key factor in understanding the 
slippage of fluids on (hydrophobic) solid surfaces. De- 
pending on the applied pressure a dewetting transition 
is found, which conversely enhances considerably the slip 
length at the surface. 

We now try to rationalize this effect at the level of 
continuum hydrodynamics, in order to capture the de- 
pendence of the slip lengths as a function of the perio- 
dicity of the roughness. Indeed, in the dewetted state 
the liquid/ wall interface takes the form of a succession of 
stripes, alternating a liquid-solid interface with a liquid- 
vapor interface. At the most simple level, the surface can 
be modelled as stripes with different slip lengths. To this 
end, we develop a macroscopic description of the hydro- 
dynamic flow on a microscopically heterogeneous surface, 
characterized by a spatially varying slip length. 



4 



III. MACROSCOPIC APPROACH : FLOW 
OVER A SURFACE WITH A HETEROGENEOUS 
MICROSCOPIC BOUNDARY CONDITION 

A. Introduction 

We study here the influence of a microscopic hetero- 
geneous h.b.c. on the macroscopic scale. Our aim is to 
define a macroscopic slip length corresponding to a he- 
terogeneous microscopic h.b.c, with a spatially varying 
value of the slip length. We therefore consider a semi- 
infinite situation, where a shear flow, with a given shear 
rate 7, is imposed at infinity. Using hydrodynamics, we 
show that far from the heterogeneous interface, that is to 
say for distances larger than the characteristic size of the 
structure, an effective slip length can be defined, which 
depends on the underlying "microscopic" slip heteroge- 
neity. 

We model the composite surface by a planar surface 
(Oxy plane) characterized by a local slip length b(x, y) 
function of the in plane coordinates (x,y). Relying on 
the microscopic results of Molecular Dynamics, one may 
assume for instance that b is locally very high over a 
vapor bubble and much lower elsewhere (for simple mo- 
dels, we may consider respectively infinite and zero slip 
lengths). The characteristics of the flow far away from 
the surface can then be determined by solving the hy- 
drodynamic equations with the h.b.c. locally given by 
b{x,y). Such an approach has been previously used to 
study local slipping heterogeneities of a simple geometric 
shape : periodic stripes parallel or perpendicular to the 
flow, corresponding to the alternance of a no slip bound- 
ary condition and an infinite slip. Philip has studied 
stripes parallel to the flow in a circular pipe of radius 
R [22]. When R — > 00, this configuration is the same as 
considering a flow at infinite over a plane, which is the 
situation studied here. In this limit, he obtains for the 
macroscopic slip length 
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where L is the periodicity of the pattern and Q is the 
fraction of the surface where the slip length is infinite. 

Lauga and Stone have established the expression of the 
macroscopic slip length for similar conditions than Philip 
(flow in a pipe of radius R with £ the fraction of infinite 
slip), but with a pattern perpendicular to the flow [23j. 
When R — > 00, their expression for the macroscopic slip 
length reduces to 
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We are interested here in a more general situation, in- 
volving an arbitrary pattern (both in geometry and ma- 
gnitude) of the microscopic slip length. 



B. Geometry and hydrodynamic equations 

We consider an incompressible simple liquid occupying 
the half-space z > 0, with an imposed shear stress 7 in 
the x direction and flowing over a fixed solid plane located 
in z — as shown in figured 

We now consider a stationary regime at low Reynolds 
number so that the Navier-Stokes equations can be writ- 
ten as 



t)V 2 U = V P 
V • (U) =0 



(7) 



where U is the velocity vector, P is the pressure and n 
the viscosity of the liquid. The boundary conditions are 
defined as follows : 

• at z = : a given local slip boundary condition 
b(x,y) which is a function of the in plane coordinate. 

U x (x;y;z = 0) 



f Hx,y) dl 
b(x,y) 



dz 1 2=0 

9Uy 

9z z=0 



= U y (x;y;z = 0) 



(8) 



t U z {x;y;z = 0) =0 



• for z 



00 



<9U 

dz 
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The effective slip length is defined in terms of the 
asymptotic velocity field, U, solution of the above equa- 
tions. In the z — » 00 limit, one expects 



U = \J S + 7ze x = (a + A fz)e x + /3e y , 



(10) 



where Us defines the "effective" velocity slip on the 
wall, with components {a, /?}: Us is a constant velocity 
vector (independent of x and y), parallel to the wall, 
which defines the macroscopic slip vector at the wall. 
The hydrodynamic equations Q and the boundary 
conditions being linear, Us must be proportional to 
7. The macroscopic slip lengths are then defined as 
B x = a/7 and B y = fl/j, here associated with a flow in 
the x direction. 

To solve this problem we introduce a new vector 
field 



V = U jze x 
and eliminate the pressure by using the vorticity 

w = V x V 
The system Q can be rewritten as : 



V 2 w 
V-(V) 









(11) 
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(13) 
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We consider a periodic pattern, with periodicity L in 
both directions. Dimcnsionless variables (starred in the 
notations) are defined as 



x — ; v 



2tt 



v 7t~ ; z 



2tt 



(14) 
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To simplify notations, we will drop the stars on the 
variables in the following, and all length and velocity 
variables must be understood as being rescaled by the 
periodicity of the pattern. The description of the reso- 
lution of equations (|13J) with h.b.c. © is given in the 
appendix, where we show that the system can be conve- 
niently rewritten in terms of the 2D Fourier components 
of the various fields in the x, y plane (we do not consider 
the Fourier transform in the z direction): 
For qjt : 



boui v 



(q) + 6(q) 



(q) 



2q5 (ffx&^q) - <? y £> x (q)) 



For q — and z = 



b o Cj, 



(q = 0) + 6(q = 0) = a 



bouj x (q = 0) 



(15) 



(16) 



where the notation "o" stands for the convolution pro- 
duct. 

The previous system for Q_ x (q) and Cj_ y (q) , in equation 
(|15fl . is solved numerically via a matrix inversion. The 
macroscopic slip lengths, B x and B y , defined in terms of 
a and /3, are then deduced by substituting this solution 
into (|16|) . Calculation details are given in the appendix. 

In the following section we present the results obtained 
for various patterns of the microscopic hydrodynamic 
boundary condition. 



C. Alternating stripes with perfect slip and no-slip 

h.b.c. 

As sketched in figure |H1 we consider a succession of 
stripes of width a characterized by a vanishing slip length 
(bi = 0), and stripes of width L — a with an infinite slip 
length (&o — > oo). The fraction of the surface with slip- 
page is C = {L — a)/L. We can consider stripes in any 
direction with respect to the flow, but we will mainly fo- 
cus on the case of stripes parallel or perpendicular to the 



flow. This system corresponds obviously to a "model" 
experimental situation corresponding to a succession of 
grooves with vapor bubbles (slipping stripes) and rough 
strips (non slipping stripes). We present below the re- 
sults for the effective macroscopic slip length B in this 
situation. 

First we study the influence of the fraction of the slip- 
ping surface, £, on the macroscopic slip length . This case 
is a benchmark for our calculation since in this specific 
situation we can compare our results with those obtained 
by Philip j2^] and Stone and Lauga [13] in the limit of 
a cylinder with an infinite radius R — > oo [2^, [2j| . As 
shown on figure [51 our approach is in excellent agree- 
ment with their analytical results. This validates the 
semi-analytical approach developed here. 

An interesting result emerging from the calculation in 
this specific geometry, is that a small percentage of non 
slipping surface (1 — £ ~ 0) is enough to decrease consi- 
derably the effective slip : even for ( ~ 98.5%, B is only 
a slightly larger than the size of the pattern of heteroge- 
neity L (here with a shear parallel to the pattern) . 



D. Alternating stripes of no-slip and partial slip 
h.b.c. 



We now consider a different situation, with a surface 
composed of alternating non slipping stripes {b\ = 0), 
and stripes with partial slip h.b.c, characterized by a 
slip length bo. The surface fraction with slip length bo is 
again defined as £. Our semi-analytical approach allows 
to consider any value of bo, and extends therefore the 
results available in the literature, which are restricted to 
the bo — oo case. We first present the results for ( = 
50%. We have plotted in figure ITUI the evolution of the 
macroscopic slip length as a function the microscopic slip 
length b in the case of a flow parallel or perpendicular 
to the strip pattern. The main observations are : 

• For low values of bo (bo < L), the macroscopic slip 
length B increases linearly with bo (of the order 
bo/ '4 for the percentage of the slipping surface con- 
sidered). 

• For high values of bo, {bo > 10 L) B asymptoti- 
cally tends to a limiting value that is a fraction of 
the periodicity (around L/10 for a flow parallel to 
the strips and L/20 for a flow perpendicular to the 
strips). 

The important point is that in the present b\ = 
case, the macroscopic slip length B is fixed by the 
smallest of the two lengths bo and L. 

The same conclusions hold for larger values of Q. 
Figure ^] shows the results obtained for the effective 
slip length in the direction parallel to the stripe, -By, for 
a surface with a fraction £ = 97% characterized by a slip 
length bo- Due to numerical limitations, we could only 
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study slip lengths bo/L < 2. The results are similar to 
those obtained at £ = 50%. 



E. Alternating stripes of infinite slip and partial 
slip h.b.c. 

We now consider a succession of stripes with an infinite 
slip length (bo — > oo) and stripes with a finite slip length 
b\ (now with a non vanishing b\). Obviously, this con- 
figuration is expected to model the situation considered 
in the first part of this work, specifically in the "super- 
hydrophobic" state. We have shown in the first part of 
this work that the interface between the liquid and the 
rough surface takes, in this state, the form of a succes- 
sion of liquid-solid stripes, characterized by a finite slip 
length (bi ~ 20cr, with a the molecular diameter), and 
liquid- vapor stripes, a priori characterized by an infinite 
slip length. We come back to this assumption in the fol- 
lowing. 

Figure IT21 shows the evolution of the macroscopic slip 
length B as a function of the microscopic slip length b\ for 
different fractions C of the slipping surface (C is defined 
here as the fraction of the surface with an infinite micro- 
scopic slip length, the rest of the surface having a finite 
slip length bi). As shown on this figure, For b\ > L, the 
effective slip length B evolves linearly with b\ for all the 
fractions £. On the other hand, one finds that for b\ = 
or b\ <C L the macroscopic slip length is fixed by the het- 
erogeneity's periodicity L, in agreement with the result of 
the section lTlI CI (for b\ = and bo — > oo). The important 
conclusion we can draw from these results is that in this 
geometry, with stripes of infinite slip length and fi- 
nite slip length, the value of the macroscopic slip 
length is determined by the largest of the finite 
lengths : bi or L. 

A simple phenomenological model can be formulated 
to explain the linearity with b\ in the case b\ > L. We 
introduce the interfacial friction coefficient A, defined by 
writing the continuity of the tangential stress a s at the 
solid-liquid interface : 

OS = r] — = XV s 

az 

where rj is the viscosity of the liquid and V s the velocity 
slip. The interfacial friction A coefficient is then related 
to the slip length b by: 

The effective friction coefficient A = can be interpreted 
as the averaged friction over the different strips, and we 
obtain accordingly the following result for the effective 
macroscopic slip length as a function of the microscopic 



which is similar to the addition rule for resistors in pa- 
rallel. 

In the case bo — > oo, we expect : 



B = 



i-C 



(18) 



1 

B 



b 



1-C 
bi 



(17) 



We have checked the linear dependance of B with bi as 
well as the dependance with £ (table HJ and we find that 
the relation is very well verified for b± > L. It is 
important to emphasize that its validity is limited to the 
case where both the slip lengths, bo and b\ are larger 
than the roughness periodicity L. Note however that in 
practice, this relationship is valid down to b\ > 0.1L (see 
figure [13 . 

For bi <C L or bi — the friction coefficient tends 
to diverge over the strip and the model does not apply. 
Physically one may attribute this failure in this limit to 
the fact that the flow is strongly perturbed (compared 
to the asymptotic form) on the poorly slipping stripes 
(bi <C L). This is not taken into account in the simple 
argument leading to Eq. I|18|) . 



F. Comparison to Molecular Dynamics results 

We are now in a position to compare the hydrodynamic 
predictions to the MD results. We focus on the super- 
hydrophobic state where the measured slip length is much 
larger than the bare slip length b\ on the smooth surface. 
In order to make the comparison with the macroscopic 
approach, it is necessary to specify the local slip lengths 
on the top of the roughness (solid-liquid interface) and 
at the liquid-vapor interface on the grooves. 

The analysis of the local profiles of velocity suggests an 
infinite slip length at the liquid- vapor interface (bo = oo), 
over a stripe of width L — a. For the slip length on the 
solid-liquid interface, on the top of the roughness, we take 
the value of the slip length b smoot h obtained on a smooth 
wall at the same pressure (b\ — b smoo th over a strip of 
width a) . We only present the results in the case of shear 
parallel to the grooves. As we have mentioned before the 
slip length of a confined liquid depends on the pressure 
while there is no pressure dependence in the macroscopic 
approach. To introduce the pressure in the macroscopic 
approach we take into account the dependence of b\ = 
bsmooth with the pressure. This dependence is measured 
independently in the simulations and depicted e.g. in 
figure 

The macroscopic slip length B is then computed u- 
sing the macroscopic approach above, as a function of 
the pressure, in the case of alternated stripes with in- 
finite slip length (bo — oo) and stripes with slip length 
b\ = b sm0 oth(P)- We compare these predictions with 
those obtained by molecular dynamic simulations with 
the same fraction £ of surface with an infinite slip length. 
In figure ITSl we compare the slip length b obtained in MD 
simulations for various pressures with the macroscopic 
slip length B computed along the lines described above. 
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The perfectly slipping surface fraction is ( = 83.3%. The 
pressure is made dimensionless using the capillary pres- 
sure. 

A few comments are in order. 

• For low reduced pressures, P/P cap < -2, there is a 
good and even quantitative agreement between the 
MD results and hydrodynamic predictions. This 
suggests that the hydrodynamic prediction, with a 
surface patterned with stripes characterized by infi- 
nite slip lengths and finite slip lengths, does capture 
the essential ingredients of heterogeneous slippage 
in the super-hydrophobic state. 

• Nevertheless we see from figure ED that this agree- 
ment deteriorates when the pressure approaches the 
capillary pressure, P cap - In particular, the depen- 
dence of the slip length on the pressure is more 
pronounced in the case of the molecular dynamics 
simulations than in the macroscopic approach. 

A possible explanation for this difference is related to 
the geometry of the liquid-vapor interface which evolves 
as the pressure reaches P cap , an ingredient which is not 
taken into account in the macroscopic approach. Figure 
n~i1 shows the evolution of the shape of the liquid- vapor 
interface, as measured in MD simulations. As can be 
seen from this figure, this interface strongly depends on 
the pressure and takes a curved shape for pressure close 
to P ca p- This suggests that the assumption of a plane 
liquid- vapor interface in the macroscopic approach in not 
appropriate at high pressure. The difference between the 
MD and hydrodynamic results indicates that the curved 
surface induces a supplementary dissipation as compared 
to the flat liquid-vapor interface, resulting in a lower ef- 
fective slip length. Indeed, at high pressure, the liquid 
penetrates in the groove so that the flow is character- 
ized by an increase of the velocity gradients compared to 
the low pressure case. These velocity gradients lead to a 
supplementary dissipation as compared to the flat inter- 
face, which could explain the difference between the two 
approaches at pressures close to the capillary pressure. 

To conclude this section, this analysis points out the 
growing importance of the shape of the interface on the 
global friction at the interface, as pressure is increased. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have considered the effect of surface 
heterogeneity on the slippage of fluid. Two approaches 
have been followed. First, we have used MD simulations 
to measure slip lengths on rough hydrophobic surfaces. 
A dewetting transition, leading to a super-hydrophobic 
state, is exhibited at small pressure, and leads to very 
large slippage of the fluid on this composite interface. 
Our results therefore confirm that mesoscopic roughness 
at the solid liquid interface can drastically modify the 



interfacial flow properties, in the same manner as it af- 
fects static wetting properties. Then on the basis of con- 
tinuum hydrodynamics, we have proposed a macroscopic 
estimate of the effective slip length associated with a sur- 
face characterized by a heterogeneous slip length pattern. 
This (semi-analytical) approach has enabled us to deter- 
mine the important factors controlling the effective slip 
length. In particular we have characterized the effective 
slip length B for a surface with a striped microscopic 
hydrodynamic boundary condition characterized by slip 
lengths bo and b\ with periodicity L. A few "rules of 
thumbs" emerge from this analysis : 

• for bi = 0, the effective slip length B is controlled 
by the smallest of the two lengths, bo or L; 

• in the case where both slip lengths bo and b\ are 
larger than the periodicity L, the effective slip 
length B can be approximated by a "inverse law" 



with £ the surface fraction with slip length bo- In 
practice, we found that this relationship holds down 
to b , bi > jqL. 

In this paper, we present results on the effective slip 
length B associated with a surface formed by alternat- 
ing stripes of different microscopic hydrodynamic bound- 
ary condition. However our results can be extended to 
more complex slipping patterns, relevant to experimen- 
tal systems. As an example, surface patterns such as 
squares or dots can be easily implemented within our 
semi-analytical approach. Note that such microscopic 
patterns can not be solved in an analytical way. 

We have then compared the results obtained within 
the two approaches, MD results and hydrodynamic des- 
cription. While the hydrodynamic approach reproduces 
quantitatively the MD results for small pressures, a dis- 
crepancy between the two approaches is exhibited at 
large pressures, close to the dewetting transition. This 
points out the importance of the shape of the liquid- 
vapor interface in the "super-hydrophobic" state on the 
dissipation at the interface. The curvature of the liquid- 
vapor interface leads to a supplementary dissipation on 
the liquid- vapor interface. 

Eventually, as an example of application of the previ- 
ous macroscopic and molecular dynamic approaches, we 
now apply these results to design a potentially highly 
slipping surface. This is motivated by potential applica- 
tions, especially in microfluidics. The motivation to de- 
sign highly slipping surfaces is linked in particular to the 
development of flat velocity profiles within microchannel, 
aiming at reducing the dispersion in the transport of dif- 
ferent species. Such a flat profile can be obtained using 
a pressure drop flow with a highly slipping channel. The 
studies we have presented show that high slippage can be 
achieved using a hydrophobic channel with "appropriate" 
roughness. We now have to design this "appropriate" 
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roughness. For microfluidic applications, let us assume 
a periodicity L = \[im. We have to determine the per- 
centage of the surface £ = which has an infinite slip 
length. We assume that the slip length of the liquid over 
the smooth hydrophobic solid is about 20 nm. Using the 
macroscopic approach, we can represent the evolution of 
the macroscopic slip length B as a function of £. Since B 
is required to be the largest, one can naively think that we 
should have the largest possible (. However, as we have 
seen with the molecular dynamic simulations, high slip- 
page requires the liquid to be in a "super-hydrophobic" 
state, so that it can "dewet" over the roughness, giving 
infinite slip length at the liquid-vapor interface. For this, 
the pressure has to be roughly lower than the capilla- 
ry pressure (as long as the grooves are deep enough). 
This capillary pressure can be expressed in terms of £ : 
P cap = ~ 27 ^ L cos6> ■ If we consider a pressure drop expe- 
riment in a microfluidic channel, the pressure drop that 
has to be applied is typically of 0.4 ■ 10 5 Pa [8j]. As re- 
presented on figure this gives us the upper limit of 
£ (C = 0.88 in our example). Therefore, if we consider 
grooves of 100 nm width and periodicity 1 /an one may 
expect slip length as large as 600 nm (figure El- This 
description does not determine the height of the grooves, 
but the main constraint for them is to be high enough, so 
that the capillary pressure does actually tune the transi- 
tion between the normal and the super-hydrophobic state 
(of course their height is also limited by the technolog- 
ical fabrication possibilities); typically here a height of 
500 nm should be appropriate. 

In summary, use of surfaces patterned at the nanome- 
ter scale, and treated to produce a 'water repellent' like 
effect, appears to be a promising way towards devices 
that would allow flow of liquids in small size channels 
with very small flow resistance. 
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Appendix: calculation details 



• For z = 



b ( x >y) 



z=0 

V z (x;y;z = 0) =0 



V y (x;y;z = 0) 



(20) 



• For z — > oo 



V = ae x + j3e y 



(21) 



We solve this problem by expanding the various fields 
in Fourier series, assuming for simplicity the dimension- 
less variables as defined in (|14|) and a pattern with pe- 
riodicity 2ir in both the x and y direction. For that we 
use the following definition for the Fourier coefficients in 
x and y directions 



V(q x ;q y ;z) = 



1 



2*) 



-7r;7r] 



r V(x;y;z)dxdy 



(22) 

where q is a two dimensional wave vector, q = q x e x +q y e y 
with discretized values q XiV = n XjV . The vorticity u> — 
V x V can then be written in Fourier space as : 



ia V - 2& 



(23) 



^ lu z = iq x V y - iq y V x 



We can deduce from equation I1K'>[1 and from the fact 
that u) must vanish as z — > oo, that : 



( u) x (q; z) = c£) 2L (q)e~« z 
u> y {q.;z) =cjy(q)e-« 2 
Lo z (q;z) =w z {q)e~ qz 



(24) 



By taking the Fourier transform of the h.b.c. I|20|l . and 
taking into account that V z (z = 0) = 0, we obtain the 
system (|2*5|) for q ^ 0. 



Here we present in some detail the resolution of the 
stationary Stokes equation with a locally imposed slip- 
ping h.b.c. We have defined in equation i|ll|) the relative 
velocity field : 

V = U - 7ze x 
which verifies the following boundary conditions : 



bow. 



(q) + S(q) = |pr (few y (q) - ^^(q)) 

b ° (q) = ^ (fe^ y (q) - <?y^(q)) + 



The Fourier transform of l|20|l gives for q = 



b o < 



(q = 0) + b(q = 0) = a 



boio x (q = 0) =f3 



^(q) 
(25) 



(26) 
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where o denotes a convolution sum. 

By solving for ui x and ui y and replacing into 
we can determine a and (3 and then the macroscopic slip 
length. Equations l|25[) only involve lv_ x , &_ y and b. Since b 
is an input of the problem, it is then possible to determine 
Uj. and Lu y . 

This resolution of Q_ x and (o_ y can not usually be done 
analytically (except for very simple patterns of b) and a 
numerical approach is required. For that, we consider a 
discretization of the elementary cell (of size (2w) 2 in the 
real space) into N 2 nodes of coordinates : 

n ( , 27r , 2 A 

\Xk', Vk) = — 7T + k x —— : — 7T + Ky — 

y ' y ; V TV y iV/ (27) 

where fc y ) G [0, N - l] 2 

The discretisation of real space implies periodicity in the 
Fourier space with period TV. We then only consider the 
wave vectors q in the interval : 



where (J7) is a vector of size 2(N 2 -1) involving only Co_ x 
and LJ y , and (B) is a vector of size 2(7V 2 -1) that only 
involves the prescribed microscopic slip length b. The 
size of the matrix (M) is 4x(7V 2 -l) 2 , which increases 
rapidly with TV. However, for most of the situation of 
interest, (M) remains a sparse matrix with many zero 
coefficients. The density of the matrix does actually 
depend on the pattern of the microscopic h.b.c. For 
a strip pattern of microscopic h.b.c, for example, only 
3 % of the elements of the matrix' elements are non zero 
(this spareness actually results from the invariance by 
translation of the pattern in the direction of the stripes). 
This allows the implementation of efficient sparse matrix 
algorithms to solve equation (|2*5|) . We typically consider 
N up to N=256 grid points both in x and y directions. 
The matrix (M) is then inverted using a conjugated 
gradient method well suited to the inversion of matrix 
of low densities. 



(q x ;q y )e[0,(N-l)} 2 

By introducing the complex vector f2(q) = & x (q) + 
juj_ y {c\) and by making a judicious combination of the 
equations (|25|) . the problem is reduced to a matrix sys- 
tem of the form : 



( M ) ( fi 



B 



(28) 
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FIG. 1: (a) Normal pressure P (in Lennard- Jones units) ver- 
sus thickness of the cell d (in Lennard- Jones units). The pe- 
riodicity is L = 18.9 a, the height of the grooves is h — 10.7 a 
and the width of the groove is L — a = 15.6 a. (b) Trans- 
verse views of the atomic configuration for different situations. 
Atoms belonging to the liquid and solid are represented by 
points and round dots, respectively. 
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FIG. 2: Measured capillary pressure P cap in the simulation 
(in Lennard- Jones units) as a function of the inverse groove 
width 1/(1/ — a). These results corresponds to grooves with 
sufficiently large depth, h, so that the measured capillary pres- 
sure becomes independent of h. 
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FIG. 3: Velocity profile in a Couette flow (black lines) : (a) 
in the normal state (high pressure) where the liquid fills the 
grooves; (b) in the superhydrophobic state (low pressure) 
where the grooves are dewetted. The upper and lower wall 
have velocity U = ±0.3 (in Lennard- Jones units). The grey 
lines are the velocity profiles obtained for a smooth surface 
(i.e. no grooves) at the same pressure. 
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FIG. 4: Evolution of the slip length b as a function of the 
distance d between the walls. The triangles (A) represent 
the values of the slip lengths on the smooth wall, the stars 
(*) represent the values of the slip lengths on the rough wall, 
characterized by periodic grooves of height h — 10.7 a, width 
a = 5 a, periodicity L = 19.7 a and cfs = 0.5 (Lennard- 
Jones units). The flow is parallel to the grooves. The circles 
(o) indicate the value of the pressure (scale on the right of the 
figure). The lines are guides for the eye. 
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FIG. 5: Evolution of the slip length (in Lennard- Jones units) 
as a function of P / PcapMary for the same geometry as in fi- 
gure 2] 
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FIG. 6: Evolution of the slip length 6 as a function of P 
(Lennard- Jones units) for grooves of different width L — a and 
same height h = 9.6 a. We only represent the slip lengths in 
the "super-hydrophobic" states. 
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FIG. 7: Schematic representation of the velocity field of the 
liquid for z > 0, for a flow in the x direction. 
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FIG. 9: We consider an succession of stripes of width a where 
the microscopic slip length is equal to zero (bi = 0) and 
strips of width L — a with an infinite microscopic slip length 
(bo — > oo) and represent the evolution of the macroscopic slip 
length B as a function of £ (£ is the fraction of the surface 
with a microscopic infinite slippage, f = (L — a)/L). The 
curves represent the analytical calculations and the dots the 
numerical resolution (top: flow perpendicular to the strips, 
bottom : flow parallel to the strips.) 
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FIG. 8: Schematic representation of the microscopic hydro- 
dynamic boundary condition, alternance of strips of no-slip 
and perfect slip. 



12 



50 
40 

n 

20 
10 



L 




• 


• • 


* 


: 


■ 










* 






] 























| 















.1 ] 10 100 



b©/L 
(a) 





0.12 




a io 




aos 


:- -. 


ao6 


m 










0.02 



I I I I I I I * III 



• * 



0.1 I 10 100 



FIG. 10: Macroscopic slip length B on a surface composed 
of alternating stripes with slip length bo and vanishing slip 
length, for £ = 50%. On figure (a), the shear flow is perpen- 
dicular to the slipping stripes; on figure (b), the shear flow 
is parallel to the slipping stripes (note the change of vertical 
scale). 
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FIG. 11: Macroscopic slip length B on a surface composed 
of alternating stripes with slip length bo and vanishing slip 
length, as a function of the microscopic slip length &o> for 
C = 97%. The shear flow is parallel the slipping strips. 
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FIG. 12: Macroscopic slip length Bona surface composed of 
alternating stripes with slip length bi and infinite slip length, 
as a function of the microscopic slip length bi . B is computed 
for a shear parallel the slipping strips and for different values 
of (. Dotted lines are a guide for the eye. In the inset the 
same points are plotted on a linear scale. The dashed lines 
are linear fits to the measured values of B. 
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FIG. 13: Comparison between MD measures of the slip length 
in the superhydrophobic state, and the macroscopic predic- 
tion of the slip length B as a function of the pressure. The 
shear is parallel the slipping stripes. The macroscopic esti- 
mate assumes alternating stripes with finite slip length 6i and 
infinite slip length; 6i is accordingly identified with b smoo th(P) 
the pressure dependent slip length on a smooth (and flat) solid 
surface, as measured independently in MD simulations. The 
percentage of the surface with an infinite slip length (macro- 
scopic simulation) is identified as the percentage of the sur- 
face with a liquid- vapor interface (in MD simulations) : here 
£ = 83.3%. Lines are just a guide for the eyes. 



FIG. 15: Evolution of the capillary pressure as a function of ( 
{Papillary = ~ 27 '" cosfl ) for 6 = 105° and 7i„ = 72 mN/m. In 



the grey zone the pressure is higher than the capillary pres- 
sure, the liquid is in the "normal" state, it is not possible 
to have high values of the slip length. The dashed lines cor- 
responds to the pressure drop for microfluidic experiments. 
The intersection between this dashed line and the capillary 
pressure's curve gives the upper value of f, here £ = 0.9. 
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FIG. 16: (a) Evolution of the macroscopic slip length B as 
a function of £, when f % of the surface has an infinite slip 
length and (1 — £) % of the surface has a slip length of 20 nm 
(the periodicity is L = 1/im). The vertical dashed corre- 
sponds to the maximum value of £ possible as explained on 
figure [T51 (b) Kind of structure that can lead to slip lengths 
up to 600 nm. 
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FIG. 14: Transverse views of the atomic configuration for 
the same geometry as the one studied on figure 1131 for (a) 
P/P 

capillary — 

0.4 and (b) P/P 

capillary =0.9. 




"Real" value of £ 


Value of £ from the 
linear fit of B (± 2%) 
as a function of bi 


97 


96.7 


83.3 


83.2 


50 


49.2 


25 


24.8 


10 


9.9 



TABLE I: Comparison of the "real" value of the fraction £ 
and the value one get with a linear fit using the equation 1181 



